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We study by computer simulations the effect of confinement on the superfluid properties of small 
two-dimensional (2D) parahydrogen clusters. For clusters of fewer than twenty molecules, the su¬ 
perfluid response in the low temperature limit is found to remain comparable in magnitude to that 
of free clusters, within a rather wide range of depth and size of the confining well. The resilience of 
the superfiuid response is attributable to the “supersolid” character of these clusters. We investi¬ 
gate the possibility of establishing a bulk 2D superfluid “cluster crystal” phase of P-H 2 , in which a 
global superfiuid response would arise from tunnelling of molecules across adjacent unit cells. The 
computed energetics suggests that for clusters of about ten molecules, such a phase may be ther¬ 
modynamically stable against the formation of the equilibrium insulating crystal, for values of the 
cluster crystal lattice constant possibly allowing tunnelling across adjacent unit cells. 

PACS numbers: 02.70.Ss,67.40.Db,67.70.+n,68.43.-h. 


I. INTRODUCTION 

The observation of the putative superfiuid phas^ of 
condensed parahydrogen (P-H 2 ) has so far been prevented 
by its strong tendency to crystallize at low temperature, 
even in reduced dimensions .12121 There exists, however, 
experimentaP evidence that hnite clusters of P-H 2 re¬ 
main liquidlike down to temperatures much lower than 
the bulk crystallization temperature, conceivably allow¬ 
ing one to probe their predicted superfiuid behaviour, 
expected to manifest itself at a tenmerature T ~ 1 K, for 
clusters of thirty molecules or less.™^ The question thus 
arises of whether one might be able to observe a macro¬ 
scopic superfiuid response in a network of interconnected 
superfiuid clusters, in which global phase coherence could 
be established by tunnelling of molecules across adjacent 
clusters. This is, in a sense, analogous to the physics of 
the recently proposed supersolid phase of soft core bosons 

in 

The above scenario could be realized experimentally, 
for example, by adsorbing P-H 2 inside a porous mate¬ 
rial such as vycor. A second avenue, perhaps affording 
greater control, exploits the predicte d sup erfiuid response 
of P-H 2 clusters confined to quasi 2D. b^ l ^^ lThe idea is that 
of fashioning a planar substrate capable of adsorbing p- 
H 2 molecules at specific sites, arranged on a regular trian¬ 
gular lattice, each designed to accommodate a number of 
molecules corresponding to a strong superfiuid response 
at low temperature (i.e., around twentjff^^. The distance 
between nearest-neighboring clusters should be chosen 
to render such a cluster crystal energetically favorable 
with respect to the formation of the ordinary uniform, 
non-superfluid crystalline phase (with just one molecule 
per unit cell), while concurrently allowing tunnelling of 
molecules across nearest neighboring wells, each one act¬ 
ing in a sense like a superfiuid quantum dot. Whether 
all of these conditions can be met is not a priori obvious. 


and furnishing a quantitative answer is the goal of this 
work. 

Besides the energetics, a second aspect to assess quan¬ 
titatively is the effect of confinement on the superfiuid 
properties of the individual clusters, which, besides be¬ 
ing obviously crucial to the goal of stabilizing the bulk 
superfiuid phase described above, is of interest in its own 
right, and might be probed experimentally, for example 
by trapping small clusters at adsorption sites of corru¬ 
gated substrates On the one hand, spatial confinement 
is expected to bring about a reduction of the superfiuid 
response of a structureless, liquidlike droplet, owing to 
the ensuing increased particle localization. However, in 
a previous studjff^ we established that 2D parahydrogen 
clusters of less than thirty molecules, turning superfiuid 
at a temperature of the order of 1 K, display a rather 
marked “supersolid” character (obviously such a defini¬ 
tion is necessarily loose, given that we are talking about 
a finite system); that is, although exchanges of identical 
molecules are frequent at low temperature, concurrently 
molecules are nonetheless spatially localized and form 
orderly structures. Thus, owing to their greater rigid¬ 
ity, superfluidity in these finite clusters may be robust 
against external confinement - more so than in helium 
droplets, for instance. Furthermore, it has been recently 
shown that confinement can actually have an enhancing 
effect on the superfiuid response of P-H 2 clusters in three 
dimensions.^ 

In order to investigate quantitatively the effect of con¬ 
finement on its structure, energetics and superfiuid prop¬ 
erties, we have carried out first principle Quantum Monte 
Carlo simulations of a single, spatially confined 2D P-H 2 
cluster at low temperature (T = 0.25 K). We studied 
clusters comprising up to thirty molecules. Confinement 
in this study is described by means of a simple gaussian 
well of varying size (typically of the order of a few A) 
and depth (up to 100 K). 
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The comparison of the computed superfluid response 
of the confined clusters with that of the corresponding 
free ones shows that, while as expected superfluidity is 
suppressed in sufficiently deep wells, nonetheless clusters 
retain their structure and superfluid response within a 
rather wide range of confinement parameters, i.e., they 
are relatively unaffected by the confinement. Suppres¬ 
sion of superfluidity takes place gradually as the well is 
deepened, affecting primarily the largest clusters. On 
the other hand, clusters with ~ 15 molecules or less re¬ 
main superfluid, at the low temperature considered here, 
even when confined in fairly deep wells. Our physical 
conclusion is that, even making allowance for the sim¬ 
plicity of the model utilized, superfluid appears to be re¬ 
markably resilient in these intriguing few-body systems; 
phrased alternatively, the quantitative requirements on 
the strength and size of the confining well may not be 
particularly stringent, at least in terms of ensuring a sig¬ 
nificant, possibly observable superfluid response of con¬ 
fined clusters. 

The energetics of the system, however, is such that, 
no matter what values of depth and size of the attrac¬ 
tive well one chooses, the distance between adjacent wells 
must be taken rather large (close to twice the size of an 
individual cluster), in order for the low-density cluster 
crystal to be energetically stable against the formation of 
the equilibrium 2D crystal. Consequently, in such a clus¬ 
ter crystal tunnelling of molecules across adjacent sites, 
necessary to establish a global superfluid response, will 
be essentially absent, for pract ical p urposes. Thus, much 
like others previously explored,this approach to sta¬ 
bilize a bulk superfluid phase of P-H 2 appears unlikely to 
succeed. 

The remainder of this article is organized as follows: in 
the next section we describe the microscopic model un¬ 
derlying the calculation and furnish basic computational 
details; we devote Sec. |ni]to a thorough illustration of 
our results, discussing the physical conclusions in Sec. 

EYl 

II. MODEL AND METHODOLOGY 

Our system of interest is modeled as a collection of N 
parahydrogen (P-H 2 ) molecules, regarded as point parti¬ 
cles of spin zero, moving in 2D in the presence of a con¬ 
fining potential. The quantum mechanical many-body 
Hamiltonian is the given by 

^ + ( 1 ) 

i=l i<j i 

where h‘^/2m = 12.031 KA^, Ti is the position of the ith 
P-H 2 molecule, = jr^ — rj|, u is the potential describ¬ 
ing the interaction of a pair of molecules, while V is the 


confining potential. We model V by means of a simple 
Gaussian well, centered at the origin, namely 

y(r) =-A exp^-^^ (2) 

While a direct experimental realization of such a con¬ 
fining potential may not be straightforward (not to our 
knowledge anyway), it contains nonetheless all the rel¬ 
evant ingredients to afford qualitative insight into the 
physics of the system, with a small number of parameters, 
thereby rendering the interpretation of the result s easie r. 
We use the well-known Silvera-Goldman potentiaP^EH to 
describe the pair-wise interaction among P-H 2 molecules. 

The thermodynamics and the structural properties of 
the above system at low temperature (r=0.25 K) have 
been studied by means of Quantum Monte Carlo simu¬ 
lations based on the continuous-space Worm Algorithm. 
Since this well-established methodology is extensively de¬ 
scribed elsewhere, we shall not review it h ere, re ferring 
instead the readers to the original references.ESESI Details 
of the calculation are standards, analogous to those em¬ 
ployed, for insta nce, i n the simulation of trapped dipo¬ 
lar Bose systems.Enm Besides the energy per molecule, 
we compute radial density profiles with respect to the 
centre of the well, as well as global and local superfluid 
response of the clusters, using standard estimators for 

finite systems 


III. RESULTS 
A. Superfluidity 

We begin by illustrating the results of our study for 
clusters trapped inside a well of size a=3 A (Eq. [^, i. e., 
roughly the radius of the inner shell of the clusters.EZl 
Figure shows the superfluid fraction Ps(W) at a tem¬ 
perature T=0.25 K, for clusters comprising up to N=30 
molecules, for wells of depth A=40, 60 and 100 K respec¬ 
tively. Also shown are the corresponding results for free 
clusters (i.e., A=0), from Ref. [T^ 

As discussed therein, a rather sharp demarcation ex¬ 
ists for free clusters, in terms of superfluid response. For, 
those with less than 26 molecules are essentially 100% 
superfluid at this low T (with the sole anomaly of A=24 
for which the superfluid response is approximately one 
half), whereas the superfluid response is suppressed in 
larger clusters. 

For a relatively shallow well {A ^ 20 K), the superfluid 
response of the confined clusters remains close (within ^ 
10%) to that of the free ones. As the depth of the con¬ 
fining well is increased, superfluidity is gradually sup¬ 
pressed, but the smallest clusters, namely those with 
N < 15, retain their superfluid properties (those with 
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FIG. 1: Superfluid fraction of 2D P-H 2 clusters confined in 
a Gaussian well of size (t=3 A and depth ^4=40 K (circles), 
A=60 K (diamonds) and A=100 K (triangles), at a tempera¬ 
ture T=0.25 K. Open triangles show results for free clusters. 
When not shown, statistical errors are at the most equal to 
symbol size. 


10 essentially entirely), even for the deepest well con¬ 
sidered here, namely with A=100 K. As shown in Fig. 
clusters whose superfluid response is most significantly 
affected by confinement are the largest ones, i.e., those 
with N > 18. 

This result may seem counterintuitive, as one might 
expect confinement to have a more disruptive effect on 
the superfluidity of smaller clusters. The reasoning would 
be that, as the well depth is increased, the molecules in 
the inner part of the cluster become localized, with the 
ensuing suppression of quantum exchanges, and thus of 
superfluidity, which might remain confined to the outer 
region of a larger cluster, where molecules enjoy greater 
mobility and where the effect of the confining potential 
is weak, for a well of size 3 A. 

This is indeed what happens, as qualitatively shown 
in Fig. [^for a cluster with N=20 molecules. Configura¬ 
tional snapshots for a free cluster (left), and one trapped 
inside a well of depth A=100 K (right), clearly show a 
much greater localization of molecules in the center of the 
cluster; this is more quantitatively illustrated by the ra¬ 
dial density profile, computed with respect to the center 
of the well, shown in Fig. Hardly any exchanges of p- 
H 2 molecules take place in the center of such a deep well. 
Exchanges still occur in the outer shell, but superfluidity 
is nonetheless suppressed to statistical noise level in the 
confined cluster, while it is nearly 100% in the free clus¬ 
ter. This is consistent with the notion that superfluidity 
in P-H 2 clusters, which have a strong “solid-like” struc¬ 
ture, crucially hinges on exchanges of molecules across 
different shells, an effect already noticed in 3D clusters 
In the presence of a confining well, superfluidity is re¬ 
silient in smaller clusters, consisting of essentially only 



FIG. 2: Configurational snapshots (particle world lines) 
yielded by a simulation of a cluster with A'^=20 P-H 2 molecules 
at r = 0.25 K. Left'. Free cluster. Right Cluster confined in¬ 
side a gaussian well of depth A = 100 K and size cr = 3 A. 
Brighter colors correspond to a higher local density. 



r(A) 


FIG. 3: Radial density profile for a cluster of A^=20 P-H 2 
molecules, confined in a Gaussian well of size (t=3 A and of 
depth j 4=100 K. Dotted line shows the corresponding profile 
for a free clusters. These profiles are computed at r=0.25 K. 
Statistical errors are not visible on the scale of the figure. 


one shell, because molecules are less compressed than in 
the case of clusters with an additional shell, and there¬ 
fore enjoy sufficient mobility, even for fairly deep wells 
(- 100 K). 

Structurally, confinement does not induce any signif¬ 
icant change, at least in the range of parameters con¬ 
sidered here. As shown in the right panel of Fig. 
molecules arrange on a 2D triangular lattice, with no 
change of structure (for example, number of molecules 
in the inner ring) for specific values of the confinement 
parameters, unlike what observed in 2D trapped dipolar 
bosons, where different solid structures are observed on 
tightening the confining well.^ 

Fig. shows ps(Af) for the same clusters as in Fig. 
but for a well of size cr=6 A. The results are quanti¬ 
tatively similar to those obtained for a tighter well, the 
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FIG. 4: Superfluid fraction of 2D P-H 2 clusters confined in 
a Gaussian well of size cr=6 A and depth A—20 K (circles), 
A=60 K (diamonds) and A=100 K (triangles), at a tempera¬ 
ture T=0.25 K. Open triangles show results for free clusters. 
When not shown, statistical errors are at the most equal to 
symbol size. 

suppression of superfluidity being more noticeable in this 
case, especially for clusters comprising between 15 and 
20 molecules, for a well of the same depth. In this case, 
the confining potential is most rapidly varying rougly be¬ 
tween the first and the second shell, for clusters of more 
than ~ 15 molecules, which has a grater suppressing ef¬ 
fect for intershell exchanges. The superfluid response of 
clusters of 15 molecules or less, on the other hand, is rel¬ 
atively unaffected for depths up to 60 K. The structure of 
the cluster is affected by the deepening of the confining 
well in the same way qualitatively illustrated in Fig. 
for a well of half the radius, i.e., molecules are increas¬ 
ingly more localized. If the characteristic radius of the 
well is further increased, essentially beyond that of the 
cluster itself, confiment becomes increasingly irrelevant, 
understandably. 

The main physical conclusion of this part of our study 
is that the superfluid response of 2D P-H 2 clusters of 
less than < 20 molecules is quantitatively rather robust 
against confinement. This is a direct consequence of 
the “supersolid” character of these clusters, which ren¬ 
ders them less compressible than liquid-like ones, con¬ 
sequently protecting their main physical properties from 
the influence of external agents. The basic physics of the 
superfluid clusters in confinement quantitatively repro¬ 
duces that of the free clusters. 


B. Energetics 

The energetics of the confined clusters is of interest in 
view of the possible stabilization of a crystal of 2D clus¬ 
ters, turning superfluid at low temperature, as explained 
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FIG. 5: Energy per P-H 2 molecule for clusters of N molecules 
(1 < < 30), trapped in a well of size ct=3 A and depth 

A=60 K (diamonds) and T^lOO K (triangles). Gircles show 
the corresponding values for free clusters (i.e., A=0). Statisti¬ 
cal errors are smaller than symbol size. Horizontal line refers 
to the ground state energy per molecule of bulk P-H 2 in its 
2D crystal equilibrium phase. 

in the Introduction. Fig. [^shows the energy per molecule 
in a cluster comprising up to thirty molecules, confined 
in a Gaussian well of varying amplitude and size cr=3 A. 
The qualitative behaviour observed in a well of twice the 
size is the same, all curves being shifted downward. 

The idea is that of “pinning” small P-H 2 clusters at 
the sites of a triangular lattice whose lattice constant d 
should be of the order of, or not much greater than, the 
characteristic size of a superfluid cluster, in order to al¬ 
low for tunnelling of outer shell molecules across adjacent 
sites. The results shown here and in Ref. [T^ suggest that 

d ~ 20 A. 

Let us assume for definiteness a number of molecules 
per unit cell N equal to 20, yielding a 2D density for 
the cluster crystal of approximately 0.058 A“^. This is 
significantly less than the CTOund state equilibrium den¬ 
sity of P-H 2 in 2D, equal to^ po = 0.0667 A“^, at which 
the system is a non-superfluid crystal with one particle 
per unit cell; the energy per molecule in such a phase is 
Co = —23.4 K. In order for the low density cluster crys¬ 
tal phase to be energetically stable against the formation 
of the ordinary 2D crystal of density po, the energy per 
P-H 2 molecule should be lower than cq + A, A being the 
average potential energy in a lattice of identical wells, of 
a given lattice constant d. This quantity can be easily 
computed numerically. 

From Fig. we see that the radius of a cluster with 20 
molecules is ^ 12 A. If the lattice constant d is taken to 
be 25 A, molecules in outer shells would have to tunnel 
across a distance of ~ I A. However, for d=25 A and (T=3 
A, we have A « —0.107 A, consistently shifting the en¬ 
ergy per particle of the 2D crystal to a lower value than 
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the energy per particle inside the corresponding well (see 
Fig. for any value of the well depth A. Thus, the 
condition of stability of the cluster crystal is not met, 
as the system finds energetically more favorable to form 
its equilibrium 2D crystal (leaving a fraction of the cell 
empty, as the density is below the equilibrium one). The 
breaking down of the cluster crystal with the formation 
of the equilibrium 2D lattice was actually observed in 
simulation. 

If the lattice constant d is taken to be nearly 30 A 
(which would entail a rather large tunnelling distance 
across sites of approximately 6 A), then A ~ —0.073 A-, 
in this case, the cluster crystal becomes energetically fa¬ 
vored for A ^ 100 K, but the superfluidity of the cluster 
is suppressed in a such a deep well, as shown in Fig. In 
order to make the cluster crystal energetically favorable, 
for a depth A such that the clusters are still superfluid, 
the lattice constant d must be taken as large as 36 A, 
making the distance across which molecules would have 
to tunnel prohibitively large. 

Increasing the width a of the well does not lead to dif¬ 
ferent physics, for a cluster of this many P-H 2 nolecules, 
because the energy per molecule in the well is shifted 
downward by an amount roughly equivalent to that by 
which the magnitude of A is increased. Moreover, the 
disruptive effect of confinement on the superfluid re¬ 
sponse is greater for this value of cr, as a result of which 
the lattice constant needed to make the cluster crystal 
thermodynamically stable is again above 30 A. 

The energy balance is more favorable for smaller clus¬ 
ters, i.e., A=10, whose radius is approximately 8 A. In 
this case, for (t=3 A, the cluster crystal with d ~ 20 A, for 
which A ^ —0.168 A, is energetically favored for A > 60 
K; it should be noted that clusters of these sizes remain 
superfluid even for such deep confining wells. Inded, on 
taking A ^ 100 K the cluster crystal is favored over the 
equilibrium 2D crystal even for d as low as ~ 19 A (be¬ 
cause we are considering molecular tunnelling, a differ¬ 
ence of 1 A is significant). Tunnelling of molecules across 
adjacent clusters would involve in this case a distance 
of 3-4 A. Whether that can allow for a superfluid phase 
at an attainable temperature remains to be established. 
The main result of this study, however, is that a super¬ 
fluid cluster crystal phase of P-H 2 , if at all attainable, 
should have a number of molecules per unit cell equal to 


ten or less. 


IV. CONCLUSIONS 

The low temperature superfluid response and ener¬ 
getics of small P-H 2 clusters trapped inside a Gaussian 
confining well have been studied by means of Quantum 
Monte Carlo simulations. The purpose was on the one 
hand to assess the robustness of the superfluid response 
predicted for the free clusters,^^ on the other that of as¬ 
sessing the possibility of stabilizing a superfluid cluster 
crystal phase of P-H 2 in 2D, analogous to that observed 
in simulations for soft core bosons.l^^ 

The main physical conclusion is that 2D clusters re¬ 
tain in confinement most of the same physical properties 
of the free systems, at least within the range of con¬ 
fining parameters explored here. Clearly, the model of 
confinement adopted here is oversimplified; a more real¬ 
istic physical model would presumably describe adsorp¬ 
tion sites as impurities around which clusters would co¬ 
alesce, i.e., with a short-distance repulsion between P-H 2 
molecules and the impurity. This may have a suppress¬ 
ing effect on the superfluid response. Also, the effect of 
foreign substitutional impurities on the superfluidity of 
the clusters has not been addressed in this study. Based 
on the computed energetics, the stabilization of a super¬ 
fluid cluster crystal phase seems possible if clusters are 
relatively small (< 10 molecules), for a lattice constant 
some 20-25% greater than the characteristic size of the 
superfluid clusters. This would require P-H 2 molexcules 
to tunnel across a distance of 3-4 A, in order for phase co¬ 
herence to be established across the whole system. Sim¬ 
ulation work aimed at ascertaining the existence of such 
a phase, as well as it possible superfluid properties, is in 
progress. 
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